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Abstract 

The principles and practice of astronomical imaging Fourier transform spectroscopy (FTS) at far-infrared wavelengths are described. The 
Mach-Zehnder interferometer design has been widely adopted for current and future imaging FTS instruments; we compare this design with 
two other common interferometer formats. Examples of three instruments based on the Mach-Zehnder design are presented. The techniques 
for retrieving astrophysical parameters from the measured spectra are discussed using calibration data obtained with the Herschel-SPIRE 
instrument. The paper concludes with an example of imaging spectroscopy obtained with the SPIRE FTS instrument. 
PACS: 95.55.-n (Astronomical and space-research instrumentation) 
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1 Introduction 

It is now well established that approximately one half of the ra- 
diant energy emitted by the universe falls in the far-infrared and 
submillimetre spectral range (30-1000 fim) IH due to two prin- 
cipal reasons. The first is that sources in the distant universe, 
ultra luminous infrared galaxies in the local universe (ULIRG), 
or protostars in our own galaxy are often shrouded in dust and 
gas. The dust efficiently scatters and absorbs shorter wavelength 
radiation from the stellar photospheres. The absorbed energy is 
subsequently re-radiated at longer wavelengths, both as contin- 
uum emission from the dust and as line emission from the ions, 
atoms and molecules present in the emitting region. In the case 
of ULIRGs, up to 99% of the radiant energy can occur in the far- 
infrared The second reason is that distant galaxies counter 
intuitively do not decrease in brightness with increasing redshift 
as (1 -H because the emission from the galaxy shifts into the 
infrared spectral range, an effect known as the negative K cor- 
rection [3 1, which compensates for cosmological dimming for 
redshifts out to z ~ 10 at far-infrared wavelengths. 

A detailed understanding of the processes underlying star for- 
mation remains as one of the outstanding questions in modern 
astrophysics. Since star formation is intimately linked with plan- 
etary formation, this has fundamental significance not only for 
understanding our immediate galactic environment, but also in 
the study of the early universe. Thus, while observations at op- 
tical wavelengths provide important information on the forma- 
tion of stars and galaxies, a detailed understanding of the phys- 
ical processes at play requires complementary observations in 
the far-infrared. Unfortunately, most of this spectral region is 
inaccessible from the ground due to absorption by the Earth's at- 
mosphere, which necessitates the use of space borne instrumen- 



tation. Moreover, since the far-infrared photon energy is < 10"^ 
eV, the instruments must be cooled, typically to liquid helium 
temperatures, in order to minimize the effects of thermal emis- 
sion from the instruments themselves, which would otherwise 
dominate the weak astronomical signal. 

In the first infrared space telescopes (Infrared Astronomical 
Satellite (IRAS) [4J, Infrared Space Observatory (ISO) |5| and 
AKARI (|6l), cooling was achieved by placing the entire tele- 
scope in a liquid helium filled cryostat, limiting the diameters 
of primary mirrors to 60 cm. While these pioneering missions 
provided our first view of the far-infrared universe, their small 
apertures resulted in relatively low spatial resolution. The Her- 
schel Space Observatory |7|, launched by the European Space 
Agency in 2009, broke this barrier by employing a 3.5 m diame- 
ter passively cooled primary mirror located outside of the instru- 
ment payload, the instrument suite being cooled to ~4 K by an 
on board supply of liquid helium. This design provides a major 
advance in spatial resolution and sensitivity, however, the sen- 
sitivity remains limited by the photon noise from the relatively 
warm (~80 K), albeit low emissivity telescope. Gains in sen- 
sitivity of two orders of magnitude will be realized by actively 
cooling larger aperture telescopes as proposed for the Japanese 
Space Agency led SPICA mission ijS). 



2 Fourier Spectroscopy 

While photometric observations in the far-infrared provide mor- 
phological information on the structure of astronomical sources, 
the composition and physical conditions of the sources can 
only be determined from spectral imaging observations. While 
different spectroscopic techniques have been explored to this 
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Figure 1; Optical configuration of the Michelson interferome- 
ter. The astronomical signal from the object plane (left) is prop- 
agated through the interferometer as described in the text, with 
the optical path difference z created by the moving mirror (right). 
The detector array is located in the image plane (bottom). 

end, imaging Fourier transform spectrometers (iFTS), have been 
widely adopted in the far-infrared. With their high throughput, 
broad spectral coverage and variable resolution, coupled with 
their well-defined instrumental line shape and intrinsic wave- 
length and intensity calibration, Fourier spectrometers have a 
long history in astronomical research Q- Furthermore, imaging 
spectroscopy is relatively easily achieved with Fourier transform 
spectrometers by the addition of an array detector, providing a 
spectrum at each spatial coordinate viewed by a pixel in the ar- 
ray. 

In the simplest type of Fourier spectrometer, the Michelson in- 
terferometer, shown in Figure[T] the incoming beam of light (e.g. 
from the telescope) is divided into two beams of equal intensity 
by a beamsplitter. After reflection from a fixed and a moving 
mirror, the beams recombine at the beamspUtter and are brought 
to a focus on the detector. 

The signal recorded by the detector as a function of the opti- 
cal path difference, z, between the recombining beams is known 
as the interferogram. In the case of a bolometer, the signal is 
proportional to the square of the electric field. The interfero- 
gram represents the autocorrelation of the electric field, Eicr), 
of the incident radiation and, in the ideal case, is given by the 
cosine Fourier transform of the spectrum, B{cr), where cr is the 
frequency expressed in wavenumbers (cm"') 
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B(cr) cos(27rzcr)dcr 



(1) 



where B{cr) is proportional to \E(o-)'^. It can be seen that in the 
ideal case the interferogram is symmetrical about the zero opti- 
cal path difference position (z - 0). In principle the spectrum 
is recovered by computing the cosine transformation of the in- 
terferogram. In practice, phase errors that arise from optical or 



electronic effects fTOl produce interferograms that are no longer 
symmetrical about z = 0, and therefore require the use of the 
complex Fourier transform. 



= J /(z)e-'2"'^Mz 



(2) 



Thus, while in principle the design of an FTS is relatively 
straightforward, involving only two interfering beams, obtaining 
the spectrum requires sophisticated mathematical analysis. 

Regardless of their design, all Fourier transform spectrome- 
ters possess two input ports and two output ports. In the case of 
a Michelson interferometer IfTTl . the input and output ports are 
superimposed. Two consequences result: the first is that 50% of 
the incident radiation is returned to the source, the second is that 
radiation emitted by the detector assembly itself can be modu- 
lated by the interferometer and return to the detector where it is 
subsequently detected. While this effect is of little or no concern 
in the optical and near-infrared, since typical instrument oper- 
ating temperatures produce no emission at these wavelengths, 
this is not the case in the far-infrared. For example, an instru- 
ment cooled to 30K has its peak emission at 100 fim. For this 
reason, Fourier transform spectrometer designs that separate the 
two input and two output ports are required for operation at long 
wavelengths. 

Historically, the Martin-Puplett (MP) polarizing interferom- 
eter 1 12 1 has been the spectrometer of choice for astronomical 
spectroscopy at far-infrared wavelengths. There are two princi- 
pal reasons: first the modulation efficiency of a polarizing beam- 
splitter is both high and uniform over a wide spectral range. Sec- 
ondly, as discussed above, the polarizing interferometer sepa- 
rates the two input and two output ports. A calibration source 
may be placed at one input port to provide absolute intensity cal- 
ibration, or the two input ports can be configured to allow differ- 
ential measurements, for example when one is trying to detect a 
weak astronomical signal in the presence of a large atmospheric 
emission component. Since the interferometric signals at the two 
output ports are complementary, they can be subtracted to double 
the modulated component of the interferogram, while in princi- 
ple cancelling any common mode noise. Alternatively, the two 
output ports can be configured to simultaneously observe two 
different wavelength ranges. One disadvantage of the MP is that, 
in its simplest form, it accepts only one polarization of the in- 
coming beam, which not only reduces the incident flux, but also 
renders the design sensitive to source polarization. Polarization 
in the continuum emission could arise from the alignment of dust 
grains in molecular clouds, or in the line emission from the pres- 
ence of magnetic fields (e.g. Zeeman effect, Goldreich-Kylafis 
effect, see 1 13 1 for a recent review). 

The Mach-Zehnder (MZ) interferometer lfT4l [Tsll . shown in 
Figure |2] exploits recent advances in the development of inten- 
sity beamsplitters ||T6|, which provide both high and uniform ef- 
ficiency over a broad spectral range. Like the MP design, the 
MZ provides access to the two input and two output ports, which 
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transform spectrometers for astronomical research at far-infrared 
wavelengths. These systems evolved from the classical Michel- 
son interferometer |fT9l to the Martin Puplett 1201 culminating in 
the Mach-Zehnder |T4l design. The following sections compare 
and contrast three MZ interferometers in which the AIG is heav- 
ily involved: The spectrometer of the SPIRE instrument on board 
the Herschel Space observatory launched in May 2009 ETIl : 
FTS-2, the spectrometer developed for use with the SCUBA-2 
camera currently being commissioned at the JCMT |22 1; and the 
spectrometer being developed for the SAFARI instrument Il23l 
on the JAXA led SPICA mission. 



Figure 2: Optical configuration of the Mach-Zehnder interferom- 
eter. After being split by the input beamsplitter, beams from both 
inputs are collimated and directed to the moving mirror assembly 
by powered mirrors. A difference in optical path length is intro- 
duced in the two beams by a translation of the moving mirror 
assembly. The collimated beams are focused by the condenser 
mirrors and recombined by the output beamsplitter; the signals at 
the two outputs represent complementary interferograms of the 
difference between the two inputs. 

can be configured in a variety of ways as discussed above. The 
principal advantage of the MZ is that all of the incident radiation 
from the astronomical source reaches the detector arrays. How- 
ever this comes at the price of increased optical path, and thus 
component size within the MZ. In order to reduce the size of the 
beamsplitters, which are typically the most demanding compo- 
nent of any interferometer, the collimating and focusing optics 
are placed within the interferometer as shown in Figure |2] 

While the separation of the two input and two output ports is 
a necessary condition to eliminate unwanted contributions to the 
measured signal from the detector environment, it is not suffi- 
cient. Each beamsplitter has a small but non-zero emissivity. It 
is thus possible for the first beamsplitter to emit light from within 
its volume, which travels through both arms of the interferome- 
ter, that results in an interference signal that is readily detectable 
with the sensitivity of modern detectors |17|. In essence, this 
is the classic double-slit experiment of quantum physics |18|. 
Since interference arising from beamsplitter emission is funda- 
mentally unavoidable, it is essential to cool the instrument to 
cryogenic temperatures. This also reduces contributions to the 
measured signal from stray light within the instrument, either 
modulated or unmodulated, which is particularly difficult to con- 
trol at long wavelengths due to the effects of diffraction. 

3 Examples of Mach-Zehnder iFTS 

For over 30 years the Astronomical Instrumentation Group 
(AIG) at the University of Lethbridge has developed Fourier 



3.1 Herschel-^VmE 

The Spectral and Photometric Imaging Receiver (SPIRE) is 
one of three scientific instruments onboard the European Space 
Agency's Herschel Space Observatory launched on 14 May 2009 
i7|. Canadian participation in SPIRE is led by the AIG. SPIRE 
contains a three-band imaging photometer operating at 250, 350 
and 500 jum, and an imaging Fourier Transform Spectrometer 
(FTS), which covers the wavelength range of 194-671 jum (447- 
1545 GHz; 14.9-51.5 cm^') |21|. The SPIRE instrument, shown 
in Figure |3] is approximately 700 x 400 x 400 mm in size and is 
supported on the ~4 K Herschel optical bench by thermally in- 
sulating mounts. The photometer and spectrometer are mounted 
on opposite sides of the SPIRE optical bench, which contains 
the optics, detector arrays (three for the photometer, and two for 
the spectrometer), an internal ^He cooler to provide the required 
detector operating temperature of ~0.3 K, filters, mechanisms, 
internal calibrators, and housekeeping thermometers. Both the 
photometer and spectrometer have cold pupil stops conjugate 
with the Herschel secondary mirror, which is the telescope sys- 
tem pupil, and defines an effective 3.29 m diameter primary aper- 
ture. All five detector arrays use hexagonally close-packed, con- 
ical feedhorn-coupled spider web neutron transmutation doped 
Germanium bolometers |f24l|. The bolometers are AC biased 
with frequency adjustable between 50 and 200 Hz, avoiding 1// 
noise from the cold JFET readout. 

The spectrometer is of the Mach-Zehnder configuration and 
uses two broadband intensity beam splitters |16| to provide spa- 
tial separation of the input and output ports. One input port views 
a 2.6' diameter field of view on the sky and the other an on- 
board reference source. Two bolometer arrays at the output ports 
cover overlapping bands of 194—313 //m (SSW) and 303-671 jum 
(SLW). As with any FTS, each scan of the moving mirror pro- 
duces an interferogram in which the spectrum of the entire band 
is encoded with the spectral resolution corresponding to the max- 
imum mirror travel. Specifications of the SPIRE spectrometer 
are summarized in Table [T] A detailed description of the SPIRE 
spectrometer is given elsewhere ETIl . 
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Table 1: Summary of the SPIRE instrument specifications. 


Snertral ranjyp 


14 Q-Si 5 cm"'- 447-l'54'5 GHz- 671-194 um 


A/taYimiim t;nf*ptral rptinliitinn 


04 cm"' - (tI \n- - 384 at fiSO iim 1 2^0 at 200 um 




T ow rO 8"^ cm"'') hitrh (0 04 cm"h 


Effective telescope diameter 


3.29 m 


JJCtliii UltliiiCLCl tlL lUUlLUIJ liiillUl 


^20 mm 


Detector optics 


Feedhom coupled 


Beam shape 


Gaussian and Multi-moded 


Angular resolution 


17-42" (FWHM equivalent beam) 


Detector arrays 


37 (SSW) + 19 (SLW) NTD Ge Bolometers 


Instantaneous field of view 


Circular, 2.6', partially vignetted 


Imaging 


2FA spacing; Nyquist sampling with 4x4 jiggle pattern 


Size 


0.7 X 0.4 X 0.4 m 


Mass 


91 kg (including three-band photometer) 


Dominant noise source 


~80 K primary telescope dish 


Operational temperature 


4.5 K (optical bench), 0.3 K (detectors), ~80 K (primary mirror) 


NEFD^ per pixel per 0.04 cm ' 


-1000 mJy (5-cr, 1 hour)'^ 


Spectral line sensitivity 


~10"'^Wm-2 (5-0-, 1 hour) 


" Noise equivalent flux density. 



A milliJansky (mJy) is a unit of flux density commonly used in astronomy; 1 Jy = 10 Wm ^ Hz 




Figure 3: A schematic of the i/eric/ieZ-SPIRE MZ Fourier trans- 
form spectrometer. Light from the Herschel telescope enters 
through input 1. The temperature of the optical bench is ~4 K; 
the short and long wavelength detector arrays (SSW, SLW) are 
at ~0.3 K. Figure adapted from Griffin et al. ||25]| . 

3.2 SCUBA-2 FTS-2 

SCUBA-2 ||26l is the first common-user instrument for sub- 
millimetre astronomy employing a large format superconduct- 
ing Transition Edge Sensor (TES) bolometer array, and is cur- 
rently operating on the James Clerk Maxwell Telescope (JCMT) 
on Mauna Kea, Hawaii. In order to achieve sky background 
hmited sensitivity under all observing conditions in the nar- 
row 450 and 850 jum atmospheric transmission windows, the 



SCUBA-2 bolometers had to be designed for low phonon noise- 
equivalent power and large saturation power, which led to the 
choice of MoCu bilayer TES devices operating at ~ 1 00 mK 1261 . 
The technological innovation of large 'CCD-like' detector arrays 
with atmospheric background noise limited performance pro- 
vides an unprecedented 2 orders of magnitude increase in map- 
ping speed over the previous camera, SCUBA fTT\. Advances in 
mapping speed and sensitivity are necessary for furthering our 
understanding of galaxy formation and evolution, a field whose 
study has seen rapid progress as a direct result of the first gener- 
ation of submillimetre instruments. Two ancillary instruments, a 
polarimeter (POL-2, 1281 ) and imaging spectrometer, have also 
been developed to further extend the capabilities of SCUBA-2. 
A Fourier transform spectrometer was selected as the optimal 
intermediate resolution spectrometer for SCUBA-2. The instru- 
ment, named FTS-2 will be primarily a galactic spectrom- 
eter (e.g. spectral index mapping of molecular clouds), but will 
also provide useful information on bright nearby galaxies and 
planetary atmospheres. FTS-2 thus fills a niche between the dual 
band SCUBA-2 continuum images and the higher spectral res- 
olution, but smaller pixel count images produced by the JCMT 
heterodyne facility instrument HARP-B OOl . 

A schematic of the FTS-2 instrument is shown in Figure |4] 
The Mach-Zehnder configuration [14] is evident, which is made 
possible by the novel intensity beam dividers. FTS-2 is config- 
ured so that both interferometer input ports are placed on the 
sky and both output ports imaged by the SCUBA-2 arrays. The 
instantaneous differencing of the two input ports cancels any 
common-mode atmospheric emission variation, which occurs on 
faster timescales than the interferogram acquisition time and is 
the dominant source of noise. Recording the complementary in- 
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Table 2: Summary of the SCUBA-2 FTS-2 instrument specifications. 



Spectral range 

Maximum spectral resolution 
Spectral resolution modes 



850 lum: 11.0-12.3 cm"'; 330-369 GHz; 910-813 fim 
450 A^m: 21.1-23.5 cm"'; 633-705 GHz; 474-426 jum 
0.006 cm-i; cr/Acr =-1950 at 850 jum, -3700 at 450 jum 
Low (0.1 cm-i), high (0.006 cm"') 



Effective telescope diameter 
Beam diameter at rooftop mirror 
Detector optics 
Beam shape 
Angular resolution 
Detector arrays 
Instantaneous field of view 
Imaging 



15 m 

-140 mm 

Direct imaging, beam defined by cold-stop 
Gaussian 

-7" FWHM at 450 yum, -14" at 850 fim 

-800 TES bolometers x 2 output ports x 2 wavebands : 

Circular, -3' at low resolution 

Q.5FA spacing at 850 jum, FA spacing at 450 fim 



3200 



Size 


-2 X 0.6 X 1.3 m 


Mass 


-600 kg 


Dominant noise source 


Atmospheric emission 


Operational temperature 


280 K (optical bench), O.IK (detectors), -270 K (primary mirror) 


NEED per pixel per 0.006 cm ' 


500 mJy at 850 jum, 2500 mJy at 450 i^m 


Spectral line sensitivity 


- 10-'** Wm-2 (5-cr, 1 hour) 




Figure 4: The FTS-2 instrument, showing the 8 aspherical fold 
mirrors, the moving back-to-back corner-cube mirror assem- 
bly, input and output beamsplitters, two retractable input port 
shutters, and retractable pickoff mirror assembly. The FTS is 
mounted perpendicular to the optical axis (which travels along 
the y axis) at the JCMT Nasmyth focus, midway through the se- 
ries of SCUBA-2 feed optics. 

terferograms in both output ports allows for gaps in the image 
due to dead pixels to be filled in. 

The telescope beam is intercepted at the intermediate Nasmyth 
image plane by a retractable pickoff mirror assembly, which di- 
verts the light through the FTS and then passes the original image 
back to the SCUBA-2 feed optics. Since the FTS was designed 
after the SCUBA-2 optics were finaUzed, compromises in the 



field of view and resolution had to be made in order to fit within 
the available space. Large aspherical diamond-turned aluminum 
mirrors are required in the interferometer arms to provide a pupil 
image at the moving mirrors, and to reimage the curved input 
image plane at the output ports. Corner-cube retroreflectors are 
necessary in the interferometer moving mirror assembly to com- 
pensate for the vertical folding of the interferometer arms. FTS-2 
is currently being commissioned at the JCMT. The key specifi- 
cations are summarized in Tablel2] A detailed description of the 
spectrometer is given elsewhere vBll . 

3.3 SPICA-SAFARI 

As noted above, despite its low emissivity, thermal emission 
from the passively cooled Herschel telescope produces a back- 
ground flux in the 20-500 jum range many orders of magnitude 
greater than that from galactic cirrus clouds, the zodiacal light 
and the 3 K cosmic microwave background. The next advance in 
far-infrared space astronomy requires a cold telescope to explore 
the cold universe. 

It is this unique low thermal background environment that the 
Japanese led Space Infrared telescope for Cosmology and Astro- 
physics (SPICA) mission will provide in the early 2020's. With 
its 3 m class primary mirror cooled to -6 K, the thermal emission 
from the SPICA telescope will be 6 orders of magnitude less than 
on JWST ll32ll or Herschel |7 |. A new generation of supercon- 
ducting transition edge sensor detectors will achieve sensitivities 
2 orders of magnitude better than existing space borne detectors 
to explore this low background environment. A consortium of 
European and Canadian scientists, which includes the authors, 
are developing a sensitive far infrared imaging spectrometer, SA- 
FARI - the SpicA FAR infrared Instrument. SAFARI will obtain 
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Figure 5: Optical design of the SPICA-SAFARI FTS. The left panel shows the Focal Plane Unit (FPU) viewed from the top; the infrared 
radiation from the telescope enters at the far left from the top as indicated by the arrow. Also shown are the two filter wheels, the optical 
beam path through the FTS (in red), and the FTS mechanism. The right panel shows the FPU from the bottom side. The light green box 
contains the input optics, the red boxes represent the Focal Plane Assemblies, the blue boxes house the filters for the readout wiring, 
and the dark green box houses the camera bay optics where the FTS output signal is directed to the three detector arrays. The middle 
panel (same orientation as the right panel) shows the optical path through the interferometer from the pick-off mirror on the bottom left 
to the three detector arrays on the top right. Figure courtesy of Roelfsema et al. Il23l . 



Table 3: Summary of the SPICA-SAFARI FTS instrument specifications. 


Spectral range 


47-294 cm-i; 1428-8817 GHz; 210-34 yum 


Maximum spectral resolution 


0.05 cm-i; cr/Acr = 2000 at 100 fim 


Spectral resolution modes 


Photometry (cr/Acr = 3), low (cr/Acr = 150-200), high (cr/Acr = 2000) 


Effective telescope diameter 


~3 m 


Beam diameter at rooftop mirror 


33 mm 


Detector optics 


Feedhom coupled 


Beam shape 


Determined by feedhom 


Angular resolution 


Diffraction limited above 40 jum (-3.5" at 40 fim, -18" at 210 yum) 


Detector arrays 


43x43 (short wave) + 34x34 (medium wave) + 18x18 (long wave) - 3329 TES detectors 


Instantaneous field of view 


Circular, ~2' 


Imaging 


0.5FA spacing at centre of each band 


Size 


TBD 


Mass 


<50kg 


Dominant noise source 


Goal zodiacal background, photon noise limited 


Operational temperature 


4.5 K (optical bench), 50 mK (detectors), 6 K (primary mirror) 


NEFD per pixel per 0.05 cm ' 


15 mJy (5-0", 1 hour) 


Spectral line sensitivity 


few 10-1'' Wm-2 (5-cr, 1 hour) 



fully spatially sampled spectra over the range 34-200 fim and 
a 2'x2' field enabling it to detect and measure, in a 900 hour 
integration, spectra of approximately one thousand times more 
extragalactic sources than the PACS instrument |33|, currently 
the most sensitive space borne spectrometer in this wavelength 
range. 

Building on the success of Herschel, the SAFARI consortium 
has selected an imaging Fourier transform spectrometer of the 



Mach-Zehnder design to achieve the spectral imaging goals. A 
schematic of the SAFARI instrument is shown in Figure |5] The 
design is described elsewhere |23|. To cover the 34-210 jum 
range with adequate spatial and spectral resolution, the inter- 
ferometer employs three detector arrays. One output port is di- 
rected to the long wavelength detector channel (1 10-210 /^m), a 
dichroic placed after the second output port splits the signal into 
the medium (60-1 10 jim) and short (34-60 jim) waveband. Full 
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Figure 6: Comparison of the measured and predicted perfor- 
mances, defined as the 5-cr spectral line sensitivity attainable in 
one hour of observing time, of the Herschel-SPIRE, SCUBA-2 
FTS-2 and SPICA-SAFARI Mach-Zehnder interferometers. For 
SAFARI, the state-of-the-art detector performance (SOAP) as- 
sumes photoconductor detectors, while the goal performance as- 
sumes TES bolometers. 1 34 1 



resolution spectra (corresponding to cr/Acr of 2000 at 100 fim) 
require a 35 mm mechanical displacement of the moving mirror 
assembly, accomplished in a time of ~2 minutes. 

The key specifications of the SAFARI instrument are summa- 
rized in Table [3] The measured and predicted sensitivities of the 
Herschel-SPIRE, SCUBA-2 FTS-2 and SPICA-SAFARI imag- 
ing Fourier transform spectrometers are compared in Figure |6] 



4 FTS Spectral Analysis 

In this section we show how a Herschel-SPIRE spectrum is an- 
alyzed to derive astrophysical parameters of the source under 
study. The analysis is generic and can be applied to any spec- 
trum obtained with a Fourier transform spectrometer Figure [7] 
shows the spectrum of one of the SPIRE calibration targets, the 
post- AGE (asymptotic giant branch) star AFGL618. This object 
is compact relative to the Herschel telescope beam and is known 
for the abundance of chemical species present in its molecular 
wind 1351 [36l : it is thus an excellent point source for calibrat- 
ing the broadband spectrometer Spacecraft pointing errors are 
potentially troublesome for a Fourier transform spectrometer as 
they can produce spectral artefacts which result in an additional 
source of error in the retrieved spectral parameters. Repeated ob- 
servations of calibration targets, such as AFGL618, throughout 
the mission, however, yield integrated line intensities that vary 
by less than 4%, which confirms that the low Herschel point- 
ing uncertainty i~2") is not a dominant source of error The 
observation shown in Figure |7] was obtained with a 10-minute 
on-source exposure [Herschel observation ID 1342242592, 2012 
March 10th]. The spectrum is the result of standard processing 



with the SPIRE FTS pipeline (v9) provided by the Herschel Sci- 
ence Centre and the SPIRE instrument control centre |.37|| . 

The spectrum reveals a baseline of continuum emission, which 
is due to dust grains. Superimposed on this continuum is an ar- 
ray of sinc-shaped spectral lines that arise from gaseous compo- 
nents in the emitting region. We have developed a line fitting 
tool for FTS spectra |38|. The algorithm performs an iterative 
fit to the entire spectrum with a superposition of a polynomial 
(for the continuum) and a number of spectrally unresolved lines. 
Since the intrinsic width of the spectral lines in AFGL618 is 
much smaller than the spectral resolution of the SPIRE FTS, the 
line shape is governed by the instrumental sine line shape, which 
has a width of 0.04 cm determined from the maximum op- 
tical path difference in the interferogram. The area (integrated 
line intensity) and peak position (central frequency) for each 
line are returned by the line fitter. By comparing the line fre- 
quencies with tabulated values of known interstellar gas species 
(JPL ||39l, CDMS EOl), lines can be identified, and are indicated 
by the labels in Figure |7] It can be seen that the rich spectrum 
of AFGL618 contains lines from molecular, atomic and ionic 
species such as CO, H2O, H, HCN, HCO+ and CCH. 

Molecular species such as CO have a ladder of rotational tran- 
sitions in the spectral range covered by SPIRE as seen in Fig- 
ure|7] Observing several transitions with varying upper level en- 
ergies, allows one to derive physical parameters for the emit- 
ting gas. A so-called 'rotation diagram' ll4Tll42ll is constructed 
under the assumptions that (a) the gas can be represented by a 
single temperature {Tx or Ty^t), (b) the lines are optically thin, 
(c) local thermodynamic equilibrium (LTE) applies, i.e., the ex- 
citation balance of the molecule is dominated by collisional ex- 
citations, and (d) the emitting region fills the telescope beam uni- 
formly. In this case, the integrated line intensity for each transi- 
tion from upper level 'u' to lower level T is directly proportional 
to A^u, the number of molecules in state 'u', which is itself a func- 
tion of Tx ' 



(rate of transitions)x(energy/photon) 



A.ihv, 



area 
1 N,„, 



4nd- 



-EJkTx 



(3) 



with h and k the Planck and Boltzmann constants, respectively, 
Aui the Einstein coefficient for spontaneous emission, Vui the fre- 
quency of the emitted photon, d the distance to the object, A^tot 
the total number of (CO) molecules summed over all states, Q 
the partition function, and the statistical weight of level u. 
Rearranging this equation yields: 



In 



guAulhVn 



-rxT,^\Qd- 



7^ ■ (4) 



A graph of y versus x is known as a 'rotation diagram'. The 
slope m of a linear fit to the data points (straight lines in Figure[8]l 
yields the gas temperature, while the intercept on the y axis (c) 
yields the total number of molecules along the line of sight. 
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Figure 7: Herschel-SPIRE spectrum of AFGL618 with identified transitions labeled. The inset diagram illustrates the fitting of two 
blended spectral lines in the observed spectrum (blue). Two individual sine profiles are shown in green and magenta, their superposition 
in red, and the residual in black, all displaced vertically for clarity. 



It can be seen from Figure [8] that the emitting region in 
AFGL618 cannot be described by a single temperature. Further, 
since the ratio of the '^CO/'^CO line fluxes varies, the '^CO lines 
are not optically thin and a full analysis requires detailed mod- 
elling of the structure of the source. Since optical depth effects 
have a greater influence on the vertical offset of the rotation di- 
agram than on its slope, it is still possible to derive meaning- 
ful rotational temperatures by invoking two temperature compo- 
nents. The temperature values indicated in Figure [8] are similar 
to those derived by Wesson et al. |36|, whose work is based on 
a SPIRE spectrum of the same source, but having a longer expo- 
sure. Even in this example of a calibration source it can be seen 
that at least two temperature components are required to fit the 
observed data. 

The rotation diagram approach outlined above can be adapted 
to correct for effects of optical depth and non-uniform beam fill- 
ing |42|. In many interstellar environments, the assumption that 
the medium is isothermal and in LTE fails and one must resort 



to a more detailed modeling of the emitting region. In general, 
such methods rely on iteration to determine interdependent pa- 
rameters: spatial distributions of dust temperature, gas temper- 
ature, molecular excitation balance and in some cases chemical 
balance ||43]| . 

The power of the Mach-Zehnder instruments described in 
this paper lies in their ability to provide astronomical imag- 
ing spectroscopy. A spectrum similar to that shown in Fig- 
ure [7] is recorded for each pixel in the detector array. To illus- 
trate the spectroscopic imaging capability of the SPIRE FTS, we 
present observations of the massive star-forming region known 
as NGC6334I ||44l|43. The Herschel observation, part of the 
Herschel Key Program "Evolution of Interstellar Dust" 1461 
[Herschel observation IDs 1342214827 (2011 Feb 26*) and 
1342251326 (2012 Sep 24*), total 6 hours], produced a spec- 
tral cube covering a 2'x2' region. Analogous to the fitting de- 
scribed above, all the CO lines were fitted in all spatial pixels 
of this cube. Figure [9] shows the integrated intensity map ob- 
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Rotation diagram for AFGL618 




"CO 



-~ J'rot = 137K 



200 300 400 

EJkiK) 




SPIRE photometer 
beam at 250 fim 



Figure 8: Rotation diagram for '^CO (blue circle points and 
solid line) and '^CO (black square points and dashed lines) in 
AFGL618. Points represent scaled measured line intensities {y 
in Equation H)) of, from left to right, the 7=4—3 up to 13-12 
transitions ('^O) and 5^ up to 14-13 ('^CO), with formal fit- 
ting uncertainties shown by vertical error bars. Straight lines 
indicate the linear regression fits to the data points. 



tained for one of these CO lines corresponding to the 7=11- 
10 transition at 42.26 cm"' (1267 GHz). The white contours 
shown in this figure indicate the distribution of dust continuum 
emission measured with the SPIRE photometer [Herschel obser- 
vation ID 1342239909], as part of the Herschel Key Program 
named HOBYS [47 1. Combining the two maps allows a direct 
comparison of the spatial extent of the (cold) dust and, in this 
case, the hot gas components of emission. Maps from all spec- 
tral lines provide astrophysicists with unique tools to compare 
and contrast the distribution of various chemical species and of 
different transitions of the same species. Combined with theoret- 
ical and numerical modeling, this leads to a better understanding 
of the physical conditions and chemical history of the regions 
under study. In addition, the determination of the central posi- 
tion of a spectral line is reliable to within a fraction (~l/20) of 
the instrumental line width (cr/Acr ai'ound 1000 for SPIRE FTS), 
yielding an uncertainty in the line centre equivalent to a velocity 
error of the order of 10 km s ' BSll. 



5 Summary 

Approximately one half of the radiant energy emitted by the uni- 
verse falls in the far-infrared/submillimetre spectral range. Sev- 
eral techniques have been developed to measure and interpret the 
information conveyed by these low energy photons. While pho- 
tometric observations provide information on the morphology of 
the source under study, the underlying physics is only revealed 



Figure 9: Integrated line intensity of CO 7=11-10 toward the 
high-mass star-forming region NGC6334I (colour scale), ob- 
served with //eric/ie/-SPIRE FTS. White contours indicate the 
distribution of dust continuum emission measured in the 250 jum 
band of SPIRE's broadband photometer. 



through spectroscopic measurements. One of the emerging spec- 
troscopic techniques, which is gaining widespread interest, is 
that of imaging Fourier transform spectroscopy (iFTS). IFTS 
provides broad spectral coverage at intermediate spectral reso- 
lution, for every pixel in the detector bolometer array. Fourier 
transformation of the measured interferograms produces a spec- 
tral data cube having two spatial and one spectral dimension. 
The designs of the classical Michelson and Martin-Puplett in- 
terferometers have been compared and contrasted to the Mach- 
Zehnder interferometer. Examples of three instruments based 
on the Mach-Zehnder design have been presented. The tech- 
niques for retrieving astrophysical parameters from measured 
FTS spectra are examined using calibration data obtained with 
the HerschelSPlKE instrument. The paper concludes with an 
example of imaging spectroscopy obtained with the SPIRE FTS 
instrument. 



Acknowledgements 

This paper is dedicated to Professor Arvid A. Schultz, friend and 
colleague, and founding member of the Physics Department at 
the University of Lethbridge. 

The authors would like to acknowledge support from Alberta 
Advanced Education and Technology, Canada Foundation for In- 
novation, Canadian Space Agency (CSA), Natural Science and 
Engineering Research Council of Canada and the University of 
Lethbridge. SPIRE has been developed by a consortium of insti- 
tutes led by Cardiff Univ. (UK) and including Univ. Lethbridge 
(Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. 



9 



Padua (Italy); lAC (Spain); Stockholm Observatory (Sweden); 
Imperial College London, RAL, UCL-MSSL, UKATC, Univ. 
Sussex (UK); Caltech, JPL, NHSC, Univ. Colorado (USA). This 
development has been supported by national funding agencies: 
CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); 
ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC, UKSA 
(UK); and NASA (USA). FTS-2 is funded through a CFI inter- 
national access award. The SAFARI consortium is led by SRON 
Netherlands Institute for Space Research, Canadian participation 
is supported by the CSA. 

References 

[1] M.G. Hauser, E. Dwek, Annu Rev Astron Astr. 39, 249 (2001) 

[2] D. Clements, L. Dunne, and S. Bales, Mon Not R Astron Soc. 403, 

274 (2010) 

[3] A.W. Blain, I. Small, R.J. Ivison, et al., Mon Not R Astron Soc. 
302, 632 (1999) 

[4] G. Ncugebauer, H.J. Habing, R. Vanduinen, et al., Astrophys J. 
278, LI (1984) 

[5] M.F. Kessler, J.A. Steinz, M.E. Anderegg, et al., Astron Astro- 
phys. 315, L27 (1996) 

[6] H. Murakami, H. Baba, R Barthel, et al., Publ Astron Soc Jpn. 59, 
S369 (2007) 

[7] G.L. Pilbratt, J.R. Riedinger, T. Passvogel, et al., Astron Astro- 
phys. 518, LI (2010) 

[8] B. Swinyard, T. Nakagawa, P. Merken, et al., Exp Astron. 23, 193 

(2009) 

[9] R Connes, Annu Rev Astron Astr. 8, 209 (1970) 

[10] S.P. Davis, M.C. Abrams, and J.W. Brault, Fourier transform spec- 
trometry. (Academic Press, 2001) 

[11] E. Hecht, Optics, 4/E. (Pearson Education, 2008) 

[12] D. Martin and E. Puplett, Infrared Physics. 10, 105 (1970) 

[13] R.M. Crutcher, Annu Rev Astron Astr. 50, 29 (2012) 

[14] P.A.R. Ade, P. A. Hamilton, and D.A. Naylor, in Fourier Transform 
Spectroscopy: New Methods and Applications, OSA Technical 

Digest. (Optical Society of America, 1999) 

[15] D.A. Naylor, B.G. Gom, I.S. Schofield, et al, Proc. SPIE Vol. 
4855, 540 (2003) 

[16] PA.R. Ade, G. Pisano, C. Tucker, et al., Proc. SPIE Vol. 6275, 

62750U (2006) 

[17] L.D. Spencer, D.A. Naylor, PA.R. Ade, et al., J Opt Soc Am A. 

28, 1805 (2011) 

[18] C. Brukner, A. Zeilinger, C. Brukner, et al., Philos. Trans. R. Soc. 
Lond. 360, 1061 (2002) 

[19] D.A. Naylor and TA. Clark, Proc. SPIE Vol. 0627, 482 (1986) 

[20] D.A. Naylor, T. Clark, and G.R. Davis, Proc. SPIE Vol. 2198, 703 
(1994) 

[21] M.J. Griffin, A. Abergel, A. Abreu, et al., Astron Astrophys. 518, 
L3 (2010) 

[22] D.A. Naylor, B.G. Gom, and B. Zhang, Proc. SPIE Vol. 6275, 

62751Z (2006) 

[23] Roelfsema, M. Giard, F. Najarro, et al., Proc. SPIE Vol. 8442 
(2012) 



[24] A.D. Turner, J.J. Bock, J.W. Beeman, et al., Appl Optics. 40, 4921 
(2001) 

[25] M. Griffin, B. Swinyard, L. Vigroux, et al., Proc. SPIE Vol. 7010, 
701006 (2008) 

[26] W.S. Holland, D. Bintley, E.L. Chapin, et al., Mon Not R Astron 

Soc. 430, 2513 (2013) 

[27] W.S. Holland, E.I. Robson, WK. Gear, et al., Mon Not R Astron 
Soc. 303, 659 (1999) 

[28] P. Bastien, E. Bissonnette, A. Simon, et al., in ASP Conf Sen 
(2011), pp. 68 

[29] D.A. Naylor and B.G. Gom, Proc. SPIE Vol. 5159, 91 (2003) 

[30] J.V. Buckle, R.E. Hills, H. Smith, et al., Mon Not R Astron Soc. 
399, 1026 (2009) 

[31] B. Gom and D. Naylor, Proc. SPIE Vol. 7741, 77412 (2010) 

[32] J.R Gardner, J.C. Mather, M. Clampin, et al.. Space Science Re- 
views. 123, 485 (2006) 

[33] A. Poglitsch, C. Waelkens, N. Geis, et al., Astron Astrophys. 518, 
L2 (2010) 

[34] D.A. Naylor, D.J. Hayton, J.V. Lindner, et al., Meas Sci Technol. 
22(2011) 

[35] R Herpin and J. Cemicharo, Astrophys J. 530, L129 (2000) 

[36] R. Wesson, J. Cemicharo, M.J. Barlow, et al., Astron Astrophys. 
518, L144 (2010) 

[37] S. Ott, in ASP Conf Ser. 434, 139 (2010) 

[38] D.A. Naylor and B.G. Gom. (in preparation) 

[39] H.M. Pickett, R.L. Poynter, E.A. Cohen, et al., J Quant Spectrosc 
Ra. 60, 883 (1998) 

[40] H.S.P Muller, R Schloder, J. Stutzki, et al., J Mol Struct. 742, 215 

(2005) 

[41] B.E. Turner, Astrophys J Suppl S. 76, 617 (1991) 

[42] PR Goldsmith and WD. Langer, Astrophys J. 517, 209 (1999) 

[43] F.F.S. van der Tak, in lAU Symposium. (201 1), pp. 449 

[44] T. Hunter, C. Brogan, S. Megeath, et al., Astrophys J. 649, 888 
(2006) 

[45] K. Qiu, F. Wyrowski, K.M. Menten, et al., Astrophys J. 743, L25 

(2011) 

[46] A. Abergel, H. Arab, M. Compigne, et al., Astron Astrophys. 518, 
L96 (2010) 

[47] F. Motte, A. Zavagno, S. Bontemps, et al., Astron Astrophys. 518, 
L77 (2010) 

[48] D.A. Naylor, J.R Baluteau, M.J. Barlow, et al., Proc. SPIE Vol. 
7731 (2010) 



10 



